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Abstract 

The mass of the W boson is determined from the direct reconstruction of 
W decays in WW — > qqqq and WW — > £uqq events in e + e~ collisions at 
LEP. The data sample corresponds to an integrated luminosity of 683 pb -1 
collected with the ALEPH detector at centre-of-mass energies up to 209 GeV. 
To minimise any effect from colour reconnection a new procedure is adopted 
in which low energy particles are not considered in the mass determination 
from the qqqq channel. The combined result from all channels is 

m w = 80.440 ±0.043(stat.)±0.024(syst.)±0.009(FSI)±0.009(LEP) GeV/c 2 , 

where FSI represents the possible effects of final state interactions in the qqqq 
channel and LEP indicates the uncertainty in the beam energy. From two- 
parameter fits to the W mass and width, the W width is found to be 

T w = 2.14 ± 0.09(stat.) ± 0.04(syst.) ± 0.05(FSI) ± 0.01(LEP) GeV. 
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1 Introduction 



The electroweak Standard Model (SM) successfully describes all interactions of quarks and 
leptons at the Z resonance provided that quantum radiative corrections are included [1]. 
In this model, the mass of the W boson (m w ) can be calculated as follows: 

mw(1 -^ ) = 71^ (1 + Ar) 

where mz (Z mass), (Fermi coupling constant) and a (fine structure constant) are 
measured with high precision. In this equation, Ar parametrises the loop corrections 
which lead to a quadratic dependence on the top quark mass, (m top ), and a weaker 
logarithmic dependence on the Higgs boson mass. A global fit of electroweak observables 
measured at the Z resonance together with the measured m top [2] yields a W mass of 
80.373 ± 0.023 GeV/c 2 [1] in the SM. 

The comparison of a direct measurement of row with this prediction was a primary 
goal of LEP, enabling a stringent test of the Standard Model to be made. This paper 
describes the final measurement of the W mass and width (r w ) from ALEPH. They are 
determined from the direct reconstruction of the invariant mass of its decay products in 
both the WW — > qqqq hadronic and WW — > 1(1 = e, fi, r)z/qq semileptonic channels. 
Measurements were published previously using the data collected at centre-of-mass (CM) 
energies of 172, 183 and 189 GeV [3, 4, 5]. The most recent ALEPH publication included 
a weighted average result obtained from the combination of all these measurements as well 
as those obtained earlier from the total W pair cross sections at 161 and 172 GeV [6, 7]. 
The statistical precision for the mass was 61 MeV/ c 2 with a systematic uncertainty of 47 
MeV/c 2 . In the last two years up to the closure of LEP in 2000, more data were collected 
at CM energies up to 209 GeV increasing the total sample by a factor of three. 

All these data, except for the small sample at 172 GeV, are included in the analysis 
corresponding to an integrated luminosity of 683 pb _1 . The data were sub-divided into 
eight samples labelled as 183, 189, 192, 196, 200, 202, 205 and 207 GeV according to their 
CM energies. This sub-division is the same as that used in the measurement of the WW 
cross section [8]. 

A constrained kinematic fit conserving energy and momentum is applied to each 
selected event in data and Monte Carlo (MC) simulation. As in previous analyses, the 
simulated mass spectra are fitted to the data using a reweighting technique to extract 
the W mass and width. Very large MC productions (> 10 6 signal events per CM energy) 
enable multi-dimensional fits to be used with significant gains in precision. The signal 
events are weighted to account for the effect of 0(a) corrections [9] in m w and r w . 

Since the statistical error on mw is now comparable with the previously published 
systematic uncertainties, a more detailed evaluation of all important uncertainties has 
been performed. In the previous analysis [5], the dominant systematic uncertainty in 
the qqqq channel was due to colour reconnection (CR). This affects the topological 
distribution of lower energy particles in an event. Two new analysis procedures have 
been adopted. In one of these, PCUT, low energy particles are not considered in the 
reconstruction of jets, whilst in the other, CONE, only particles close to the jet axes are 
used. These significantly reduce the difficult-to-estimate CR uncertainty at the expense 
of some statistical power in this channel. The effect of these modified reconstructions 
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has been checked using di-jets in the £vqq channels, where no final state interactions are 
present between the W's. 

The paper is organised as follows. In Section 2, the important properties of the 
ALEPH detector are described. In Section 3, the event reconstruction procedures and 
calibrations are recalled and detailed studies of the detector simulation reported. Section 4 
contains a full description of the event samples generated for the signal and background 
processes involved. Section 5 describes the event selection and kinematic reconstruction 
procedures in the different channels highlighting, where appropriate, the modifications 
and improvements applied since the earlier analyses at 183 and 189 GeV [4, 5]. Section 6 
describes the extraction of mw and IV- Section 7 describes the specific studies made to 
set a limit on colour reconnection from the data using event reconstructions where low 
momentum particles or particles between jets are excluded. Section 8 describes all studies 
of systematic uncertainties. The analysis of radiative returns to the Z peak is reported in 
Section 9, providing a cross check on the W mass measurement. The measurements in all 
channels are combined in Section 10, taking into account common sources of systematic 
uncertainties. The W masses measured in the purely hadronic and combined semileptonic 
channels are compared in this section. Final conclusions are given in Section 11. 



2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in Ref. [10] and of its 
performance in Ref. [11]. The tracking detectors include a silicon vertex detector 
(VDET), a cylindrical drift chamber and a large time projection chamber (TPC) which 
measures up to 31 coordinates along the charged particle trajectories. A 1.5 T axial 
magnetic field, provided by a superconducting solenoidal coil, yields a resolution of 
5pt/pt = 6 x 10 -4 pr © 0.005 (pr in GeV/c). Charged particle tracks reconstructed 
with at least four hits in the TPC and originating from within a cylinder of 2 cm radius 
and 20 cm length, centred on the nominal interaction point and parallel to the beam axis, 
are called good tracks. In addition to its role as a tracking device, the TPC also measures 
the specific energy loss by ionisation dE/dx. 

Electrons and photons are identified in the electromagnetic calorimeter (ECAL) by 
their characteristic longitudinal and transverse shower development. The calorimeter is a 
lead/ wire-plane sampling detector with fine readout segmentation. Each tower element is 
projective, subtending an angle of ~ 1° in both 9 and <p, and segmented longitudinally into 
three stacks. It provides a relative energy resolution of 0.180/\^E + 0.009 (E in GeV) 
for isolated electrons and photons. The three-dimensional fine segmentation allows a 
good spatial resolution to be achieved for photons and 7r°'s in jets. Such deposits are 
separately identified and their energies evaluated by a fine clustering algorithm [11]. 
Muons are identified by their penetration pattern in the hadron calorimeter (HCAL), 
a 1.2 m thick iron yoke instrumented with 23 layers of streamer tubes, together with two 
surrounding double layers of muon chambers. The hadron calorimeter also provides a 
measurement of the energies of charged and neutral hadrons with a relative resolution of 
0.85/ \/E (E in GeV). At low polar angles, electromagnetic energy deposits are detected 
in the luminosity calorimeters (LCAL and SiCAL [11]) down to 34 mrad with respect to 
the beam axis. 
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3 Event Reconstruction 



3.1 Simulation of electromagnetic deposits 

Aiming for a precise measurement of mw to 1 part in 2000 imposes a level of understanding 
of this detector and its simulation not required in previous analyses of ALEPH data. The 
fine granularity and longitudinal segmentation of the ECAL detector elements [10] allow 
nearby energy deposits to be identified. The treatment of these deposits has been revised 
following a detailed simulation. 

The normal simulation of the lateral and longitudinal development of electromagnetic 
energy deposits in the ECAL tower elements is based upon a parametrisation of showers 
measured in a test beam. This parametrisation was employed in the generation of all 
reference events used in the analysis (~ 10 8 events). It provides a good description 
of the individual shower cores but fails to simulate the correlated fluctuations in their 
development through the sampling layers of ECAL, which can lead to the production 
of objects separate from the main deposit, called satellites. Mostly below 1 GeV and 
confined to one stack, the observed rate of such objects significantly exceeds expectation. 
To understand the origin of this discrepancy, a more complete simulation (FULLSIM) of 
the response of ECAL to electromagnetic showers was developed using GEANT3 [12]. The 
effect of correlated fluctuations is included. As expected, there is better agreement in 
the reproduction of low energy satellites. However, since FULLSIM was restricted to an 
average medium for the ECAL sampling layers, it does not describe the lateral shower 
shape as well as the parametrisation. Consequently, its use was confined to the study of 
calorimeter systematic effects for which samples of ~ 10 6 events were generated. 

Specific studies with 45 GeV Bhabha electrons show an excess in the data of objects 
formed entirely from connected elements from within the same stack. Similar effects are 
seen in the close neighbourhood of particles in jets. Not identified as electromagnetic, 
all single stack objects are removed from both data and simulated events unless related 
to a good track or a HCAL energy deposit. After this ECAL 'cleaning' process, the 
multiplicity of single stack objects in ECAL matches the prediction from FULLSIM. The 
multiplicity of identified photons is unaffected by this procedure. The ECAL cleaning 
removes ~3% from the total energy of a hadronic jet both in the data and FULLSIM. 

3.2 Energy flow 

The total visible energy and momentum per event and thus the missing energy and 
momentum, are evaluated by an energy flow reconstruction algorithm [11] which combines 
all measurements from calorimeters and tracking devices. The algorithm also provides a 
list of charged and neutral reconstructed particles, called energy flow objects, from which 
jets are reconstructed. The four-momentum of a jet is defined as the sum of the four- 
momenta of all particles in the jet, assuming the pion mass for all charged hadrons. 
The typical jet angular resolution is 30 mrad in space. The jet energy resolution is 
approximately a E . ct = (0.Q^E~ t + 0.6) GeV where E- ]Ct (in GeV) is the jet energy. 

In order to bring better agreement between data and simulation, all energy flow objects 
in data and simulated events found to subtend angles less than 15° to the beams are 
rejected. All photonic and hadronic objects identified only in the ECAL are rejected if 
their energies are less than 1.5 GeV. Hadronic objects identified in the third stack of the 



3 



ECAL combined spatially with an HCAL deposit are rejected if their energies are less 
than 2 GeV. Objects with energies below these thresholds are not perfectly described by 
the simulation of the detectors. 

3.3 Calibrations 

Large samples of Z decays were collected at 91.2 GeV CM energy at the start and end of 
LEP2 running each year. Di-lepton and di-jet events were used to monitor the performance 
of the detector and to compare reconstructed particle and jet four-momenta with the 
predictions of the simulation. The following subsections describe the corrections applied 
where significant discrepancies between data and simulation were found. 

3.3.1 Charged particles 

For charged particles, small sagitta corrections are applied in data as determined using 
di-muons. They are proportional to momentum and opposite in sign for positively and 
negatively charged particles reaching a relative difference of 2% for 45.6 GeV/ c tracks at 
the smallest polar angles. 

3.3.2 Momentum of isolated leptons in e,/xi/qq events 

Electron candidates from semileptonic W decays are corrected for energy losses due to 
bremsstrahlung in the detector material by combining their four-momenta with those 
of any detected photons that are consistent with this hypothesis. These photons can 
appear either as an excess of energy in the ECAL electron cluster or as a separate deposit 
topologically consistent with bremsstrahlung. This correction is not applied when the 
electron is accompanied by other charged particles with summed momenta greater than 5 
GeV/ c within 6° of the electron track. In addition, for muon and electron candidates, a 
search is made for isolated final state (FSR) photons associated with the lepton. Such a 
photon must be closer to the good lepton track than to any other object or the beam axis 
and at least 40° away from any other good track. Their four-momenta are then combined. 

In addition to the treatment of sagitta distortion described in subsection 3.3.1, the 
simulation of electrons from Bhabha events at 91.2, 130-136 and 183-209 GeV was 
compared with data. Small systematic biases as a function of polar angle and electron 
energy E c were found arising from an imperfect simulation of saturation and leakage 
effects in ECAL. The main effect is a global relative shift of 0.45% parametrised as 
AE C /E C (%) = 0.45(0.04) - 4(8) • 10~ 4 [£ c - 45.6 GeV] (errors in brackets). This is 
applied as a correction to the simulated electrons to match the data. A similar study 
for muons yields a small miscalibration of the momentum at 45 GeV (0.08±0.03%) with 
no significant dependence on momentum or polar angle. In this case, no corrections are 
applied. 

3.3.3 Identified photons 

Using the energy flow algorithm, photons are identified in ECAL both in isolation and 
from within clusters of overlapping objects. Any bias in the photon energies from the 
simulated events relative to data was investigated by comparing n° mass distributions 
made from tau pairs at the Z. In addition, directly measured single photons from ///ry 
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events were compared event-by-event with the corresponding kinematically reconstructed 
values. Small biases are corrected to match the data, parametrised separately for the 
barrel, endcaps and the 'overlap' region in between. 



3.3.4 Jets 

Following these corrections to the charged hadrons and photons within jets, simulated 
hadronic events at the Z with energies of 45 GeV were compared with the data. The 
hadrons are clustered into two jets using the DURHAM-PE algorithm [13]. Only qq events 
with thrust values in the range 0.8 to 0.9925 are used, to suppress three-jet configurations 
and tau pairs. Using all Z calibration data collected during the LEP2 data taking periods, 
a statistical precision of about 0.2% on jet energies is obtained. Figure 1 shows the ratio 
of jet energies in data to simulation, determined from the mean values in each bin, as a 
function of jet polar angle cos 6^. The relative biases in the barrel region do not exceed 
0.5% and reach a maximum of 2.5% for | cos#j e t| > 0.95. The Monte Carlo reconstructed 
jet energies are corrected bin-by-bin for these biases as a function of cos 0j e t before event 
kinematic fits are applied. Figure 1 also shows the relative jet energy resolutions as 
determined from the RMS values of the distributions in each cos 0j et bin. The simulation 
agrees with the data to within 1% for the barrel and 4% for the endcaps; no correction is 
applied. 
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Figure 1: The ratios of jet energies (circles) and resolutions (triangles) at the Z peak between data 
taken at the Z and corresponding simulation as a function of cos#j Ct . 



Variations in the corrections as a function of jet energy have been studied above 45 
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GeV by comparing di-jet event samples from data and simulation at CM energies from 
130 to 209 GeV. Below 45 GeV, a large sample of jets with energies centred at 30±7 
GeV were obtained from three-jet events at the Z peak. In this way, the full range of 
jet energies from WW decays is covered. Figure 2 shows the average ratios of measured 
jet energies in data to simulation for the barrel and endcap regions separately for three 
values of jet energies from 30 to 98 GeV. No significant deviations are observed from the 
ratios found for 45 GeV jets. 
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Figure 2: Data/MC ratio of average E- ]ot /E h 
di-jet events at the Z and (c) high energy di-jet events 
| cosset | < 0.7. 



,. mn for £j et energies of (a) 3-jet events at the Z, (b) 



Closed points: |cos#j C t| > 0.7, open points: 



4 Simulation of Physics Processes 

The K0RALW event generator, version 1.51 [14], is used to produce W pair events. These 
events are weighted by the 0(a) correction to the doubly resonant W-pair process using 
the YFSWW3 program version 1.16 [9]. Within K0RALW all four-fermion (4f) diagrams 
producing WW-like final states are computed, including Cabbibo suppressed decay modes, 
using the fixed-width scheme for W and Z propagators. The JETSET 7.4 [15] or PYTHIA 
6 . 1 [16] packages are used for the hadronisation of quarks in the final states. Their 
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parameters are tuned at the Z from a selection of qq events with anti-b tagging. Colour 
reconnection and Bose-Einstein final state interactions are not included. A sample of 10 6 
4f events was generated with KORALW at each of eight CM energies ranging from 182.7 to 
206.5 GeV [8]. The W mass was set to 80.35 GeV/ c 2 and the width taken from Standard 
Model (SM) predictions to be 2.094 GeV. These samples are used as reference samples 
for fitting to the data in the reweighting procedure (Sec. 6), as well as for the study 
of detector systematic errors. Additional samples of 200k events were generated with 
W masses up to 0.5 GeV/ c 2 and W widths up to 0.6 GeV different from the reference 
sample, for checking the stability of the results. Also, an independent sample of 500k 
W pair events was generated at each CM energy with KORALW restricted to the doubly 
resonant CC03 diagrams [8]. This sample is used to train neural networks and parametrise 
the corrections used in the kinematic fitting. 

For studies of the systematic errors from fragmentation in W decays, 10 6 W pair 
events generated with KORALW were hadronised using JETSET, HERWIG 6.2 [17] and 
ARIADNE 4 . 10 [18] and processed through the detector simulation. To suppress statistical 
fluctuations in the comparison between these models, the same events at the parton level 
were used. Similarly, simulated samples of 100k to 500k events, generated with KORALW, 
were hadronised with modified versions of JETSET [19, 20], HERWIG and ARIADNE [21] 
containing various implementations of colour reconnection, to assess the influence of final 
state interactions between W decay products on the mass and width. The KORALW 
events were also hadronised with a version of JETSET that includes Bose Einstein 
correlations [22, 23], to determine their influence on the W mass and width measurements. 

Simulated samples of events of at least hundred times the data luminosity were 
generated for all background processes at each CM energy. The e + e~ — > qq(7) events 
were generated with KK version 4.14 [24] with hadronisation performed by PYTHIA and 
including final state photon radiation (FSR) in the parton shower step. Interference 
between initial and final state was not taken into account. Events from ZZ-like final 
states were generated using PYTHIA (NC08 diagrams), but particular care was taken to 
avoid double counting of ZZ events already included in the signal generation as WW-like 
events (i.e. uudd, fi + n~i>v,..). The same applies to Zee final states, generated with a 12 
GeV/ c 2 minimum mass for the Z system, for what concerns double counting of e + e~uu 
events. Two-photon (77) reactions into leptons and hadrons were simulated with the 
PH0T02 [25], PYTHIA and HERWIG generators but no events survived the selection cuts in 
the qqqq and tvqq channels. Di-lepton final states were simulated using KK for tt{j) and 
/i/i(7) and BHWIDE 1.01 [26] for ee(7) events. Table 1 lists the number of simulated events 
passed through GEANT, generated for each process at each CM energy with corresponding 
integrated luminosities for the data. 

5 Event selections and kinematic reconstruction 

In the following subsections, the event selections and kinematic reconstruction procedures 
for the mass extraction are described for the following four classes of WW events qqqq, 
ez/qq, /zz/qq and rz/qq. The selections are those required for the WW cross section 
measurement [8]. For the qqqq, ez/qq, and /iz/qq channels, the cuts developed earlier 
at 189 GeV [5] for the leptons and jets are used followed by re-optimised neural networks 
for the higher CM energies. A new selection has been developed for the rz/qq channel. 
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Table 1: Overview of the numbers of simulated events generated (in units of 1000 events) for each 
process type at eight average centre-of-mass energies and corresponding data integrated luminosities. 
Events generated from the same 4f (2f) samples are shown in italics. 



year 


1997 


1998 


1999 


1999 


1999 


1999 


2000 


2000 


Energy (GeV) 


182.65 


188.63 


191.58 


195.52 


199.52 


201.62 


204.86 


206.53 


Luminosity and 


56.81 


174.21 


28.93 


79.86 


86.28 


41.89 


81.41 


133.21 


total error (ptT 1 ) 


±0.31 


±0.77 


±0.14 


±0.36 


±0.38 


±0.20 


±0.38 


±0.60 


4/- signal 


1200 


1400 


1200 


1200 


1200 


1200 


1200 


1400 


4/-JETSET 




1000 












1000 


4/-HERWIG 




1000 












1000 


4/-ARIADNE 




1000 












1000 


ZZ (NC08) 


200 


200 


200 


200 


200 


200 


200 


200 


Zee (> 12 GeV/c 2 ) 


200 


200 


200 


200 


200 


200 


200 


200 


e + e~ 


1000 


3000 


1000 


3000 


3000 


1000 


3000 


3200 




300 


300 


300 


350 


300 


300 


300 


300 


T + T - 


100 


160 


100 


100 


100 


100 


100 


100 


qq 


2000 


1200 


2000 


2000 


2000 


2000 


2000 


1000 


qq-JETSET 




1000 












1000 


qq-HERWIG 




1000 












1000 


qq-ARIADNE 




750 












750 


qq-BE32 


150 


150 


150 


150 


150 


150 


150 


150 


77 — > e + e~ 


700 


1200 


700 


900 


900 


1100 


900 


1100 


77 -> n+fj,- 


600 


1200 


700 


900 


900 


800 


900 


1100 


77 — > t + t~ 


800 


1200 


700 


900 


900 


800 


900 


1100 


77 — > hadrons 


















un-tagged 


1000 


3000 


500 


1500 


1500 


500 


1500 


2500 


tagged 


500 


1000 


500 


1000 


1000 


500 


1000 


1000 


CC03-JETSET 


500 


500 


500 


500 


500 


500 


500 


500 


CR model SKI 


500 


500 


500 


500 


500 


500 


500 


500 


CR model SKII 




500 












500 


CR model SKII' 




500 












500 


CR model GAL 


250 


250 


250 


250 


250 


250 


250 


250 


BE32 2 models 


300 


300 


200 


200 


300 


200 


200 


300 


CC03- ARIADNE 


500 


500 


500 


500 


500 


500 


500 


500 


CR model AR2 


500 


500 


500 


500 


500 


500 


500 


500 


CC03-HERWIG 


500 


500 


500 


500 


500 


500 


500 


500 


CR model 11% 


500 


500 


500 


500 


500 


500 


500 


500 



All selections are mutually exclusive. 
5.1 WW — > qqqq selection 

A first preselection step aims at removing events with an energetic undetected initial state 
radiation (ISR) photon from radiative returns to the Z by requiring that the absolute value 
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of the total longitudinal momentum be less than 1.5(M vis — Mz) where M vis is the observed 
visible mass. All accepted particles are then forced to form four jets using the DURHAM-PE 
algorithm [13]. Only events where the jet resolution parameter, 1/34, is larger than 0.001 
are kept. To reject qq events with a visible ISR photon, none of the four jets can have 
more than 95% of electromagnetic energy in a 1° cone around any particle included in the 
jet. Four-fermion final states in which one of the fermions is a charged lepton are rejected 
by requiring that the leading charged particle of each jet carries less than 90% of the jet 
energy. 

The same neural network (NN) as in Ref. [8], trained at five CM energies (189, 196, 
200, 205 and 207 GeV) on Monte Carlo CC03 and background events, is used to tag the 
preselected events. There are 14 input variables based on global event properties, heavy 
quark flavour tagging, reconstructed jet properties and WW kinematics. The signal is well 
separated from the qq(7) background with 90% efficiency and 80% purity by requiring a 
NN output in excess of 0.3 [8]. 

According to the simulation, a significant fraction (~6%) of the accepted events are 
accompanied by an ISR photon that can be detected in the calorimeters separately from 
the hadronic jets. Such photons can be removed from the jet clustering process, thus 
improving the invariant mass resolution for W pairs. Studies at 189 GeV show that such 
photons with energies above 3 GeV are identified in SiCAL or LCAL and above 5 GeV 
in ECAL with an overall efficiency of 63% and purity of 72% if an isolation criterion 
based on a minimum angular separation from the closest energy flow object is applied. 
The minimum separation applied is 8° in SiCAL or LCAL and 18° in ECAL for all CM 
energies. These events are treated differently in the subsequent kinematic fit. 

5.2 WW — > ivqq selection 

A preselection common to the three lepton topologies requires at least seven tracks in the 
event. Background from qq events is reduced by requiring the estimated sum of missing 
energy and missing momentum to be greater than 35 GeV. The Z7 events in which the 
photon is undetected are rejected by requiring the missing longitudinal momentum to be 
smaller than 

Max((s - M|) /2v^ - 27.5 GeV, (^/s - M 2 z /^s - \j - pV - 6 GeV) 

where is the transverse missing momentum and ft is the missing energy. 

Following the identification of the lepton and associated objects, the remaining 
particles are clustered into two jets using the DURHAM-PE algorithm as in the qqqq channel. 

5.2.1 ei/qq and /^i/qq selection 

In addition to the common preselection, a tighter cut is used on the total visible energy 
and visible longitudinal momentum to further reject Z7 events: 

E**( s - Ml)/{s + Ml) - > 5 GeV 

where E V1S and PJ 1S are the visible energy and longitudinal momentum, respectively. 

The lepton candidate is chosen as the good track with the largest P sin (fly /2) where 
P is the track momentum and #ij is the angle from the track to the closest jet clustered 
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from the remaining tracks using the DURHAM-PE algorithm (y cut = 0.0003). Events are 
further considered if this lepton candidate satisfies the electron or muon criteria defined 
in Ref. [8] and if the sum of the lepton and missing energies is greater than 30 GeV. 

Two different NN's have been trained to select and classify ez/qq and /zz/qq signal 
events [8]. Both use three discriminant variables, the event transverse momentum, the 
lepton energy and the lepton isolation. The event is classified as ez/qq or /iz/qq if the 
corresponding NN output value is larger than 0.6 [8]. The efficiency and purity of the 
ez/qq selection are 82% and 93% respectively. The corresponding values for the /xz/qq 
channel are 89% and 98%. 

Detailed studies of neutral objects not already classified as bremsstrahlung within 2.5° 
of the electron track impact point on ECAL show a higher multiplicity than expected 
even after the removal of single stack objects (Sec. 3.1). The reference simulation fails 
to reproduce the data for angles up to 8°. Further studies show that a smaller but still 
significant excess of charged objects are present in the data for both ez/qq and /zz/qq 
events. Although the summed energy of these objects near the isolated lepton is small, 
their impact on the closest jet is significant, especially for the ez/qq channel. Therefore, 
all these objects up to 8° from the lepton are removed from the jet reconstruction. Also, 
they are not included in the calculation of the lepton four-momentum. 

5.2.2 ri/qq selection 

A new selection has been designed [8], based on an improved tau reconstruction [27]. 
Leptonic tau decays are searched for by examining those events with e or /i candidates 
which fail the ez/qq or /iz/qq selection. These events are subjected to a similar three 
variable NN but trained on leptonic tau decays. Events with the NN output greater than 
0.4 are kept [8]. 

After removing the events which have satisfied any of the three variable NN selections 
for ez/qq, /iz/qq or rz/qq, the remaining events are further examined for additional rz/qq 
final states. Use is made of the fact that one-prong tau decays are characterised by a low 
visible mass with a mean about 0.75 GeV/ c 2 . The first step is to perform a jet clustering 
using the JADE algorithm [28] with a low y cut = (0.75/£ vis ) 2 (E vis in GeV). The tau 
candidate is defined as the jet which maximises pj (1 — cos#j), where 9j is the smallest 
angle with respect to other jets and pj is the jet momentum. The event is then subjected 
to additional cuts, in particular the invariant mass of the hadronic recoil system to the 
tau candidate must be in the range 60 to 105 GeV/c 2 . For those events which fail, the 
procedure is repeated with increasingly higher values of y cut . When this exceeds (5.0/E vis ) 2 
the iterations are stopped and the event is kept requiring only that the recoiling invariant 
mass is larger than 20 GeV/ c 2 [8]. 

If a r-jet candidate is found, the event is subjected to further cuts to remove the main 
backgrounds. Most of the 77 interactions are rejected by requiring the visible mass of 
the event to be larger than 50 GeV/c 2 and the missing transverse momentum greater 
than 10 GeV/c. The event is divided into two hemispheres with respect to a plane 
perpendicular to the thrust axis. The acollinearity angle between the two hemispheres 
is required to be less than 175° to reject most of the qq background. About 80% of the 
events with a tau candidate satisfy these cuts but significant background remains, mainly 
from qq events. These events are then subjected to a 15 variable neural network. The 
event is selected if the result is greater than 0.4. The efficiency and purity of the rz/qq 
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selection are 65% and 86% respectively [8]. 



5.3 Kinematic fit 

The biases and resolutions used in the kinematic fits for the jet energies and directions 
are determined from an independent CC03 simulated sample. The distributions of the 
differences between the reconstructed jet energies and angles and those of the jets built 
directly from the generated particles are binned in jet energy and polar angle. Each of 
these distributions is fitted to a Gaussian and the mean values and sigmas are fed to the 
fitting algorithms. 

Except for the rz/qq channel, W pair events are treated as four body final states with 
either four jets or two jets, a charged lepton and neutrino to which the measured missing 
momentum is assigned. For each selected event, two invariant masses are computed from 
the W decay products. In order to improve resolution, kinematic fits are made with 
the constraint of event four momentum conservation and fixing the velocities (p/E) of 
the jets to their measured values. Imposing energy and momentum conservation alone 
corresponds to a four-constraint (4C) fit in the case of fully hadronic events, and a one- 
constraint (1C) fit in the case of semileptonic events, giving two different fitted masses 
per event. An equal mass constraint for the two bosons corresponds, respectively, to a 
five (5C) or two-constraint (2C) fit. In the rz/qq channel, since the tau energy is largely 
unknown due to neutrinos in the tau decay, only the hadronic side of the event is used 
with the sole constraint of the beam energy. 

The average raw resolution of 12% on the total jet momentum improves by a factor 2 
and by a factor up to 5 for polar angles down to 20 degrees, due largely to the kinematic 
fitting. 

For all classes of events the fits converge successfully producing flat x 2 probability 
distributions for P(x 2 ) > 0.05. The peak below P(x 2 ) — 0.05 is populated by events that 
do not fully satisfy the fitting hypothesis. Monte Carlo studies show that approximately 
half of these events have ISR energies greater than 0.5 GeV, leading to a significant 
positive bias in the reconstructed di-jet masses. However, these events are not removed 
since the simulation adequately describes the observed x 2 probability distributions in all 
channels. 

In the qqqq channel for those events with an identified ISR photon in the detector, 
the procedure of event clustering and fitting is modified [5]. In this case, the energy flow 
objects from which the ISR photon has been removed are forced into four jets. The fit is 
performed taking into account the modified constraints 



E(i*»fl) = (vM 



J2(EuPi) = {\fs-E 1 ,-p 1 ) 



i=l 



Of the 4861 data events selected after all cuts, 220 are treated in this way. Monte Carlo 
studies at 189 GeV show that the invariant mass resolution for these events improves 
from 4.1 to 2.9 GeV/c 2 and the mean displacement of the masses from their true values 
is zero within error. The improvement in the expected error on mw for all selected events 
is -2%. 
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Jet pairing in the qqqq channel 

At most one of the three possible jet pairings is chosen, based on the the CC03 matrix 
element \M.{pi 1 ,P{ 2 ,p^,P{ ll , m w )| 2 , where the p^s denote the kinematically fitted four- 
momenta of the respective jets and m w f the reference W mass, taken to be 80.35 GeV/ c 2 . 
The combination with the largest value of \M\ 2 is in general selected (in 90% of the cases), 
provided that (a) it does not have the smallest sum of jet-jet angles and (b) both fitted 
masses lie in the [60,110] GeV/c 2 window. Otherwise (in 10% of the cases) if it satisfies 
the same criteria, the combination with the next-to-largest value of \Ai\ 2 is chosen. If the 
pairings with the two largest values of \Ai\ 2 are not accepted, the event is rejected. At 
189 GeV [5] for example, the fraction of kinematically fitted signal events surviving these 
criteria is 80%. Of these events, 90% are found to have the correct combination of di-jets 
when comparing their directions to those of the original W — > qq decays. The bias from 
the choice of reference mass is found to be negligible. In addition, the combinatorial and 
physical backgrounds do not show particular structure in the defined mass window. 

6 Extraction of the W mass and width 

The W boson mass and width are extracted by fitting simulated invariant mass spectra 
to the observed distributions. As in previous analyses [3, 4, 5] an unbinned maximum 
likelihood procedure is employed to find the best fits, using probability density functions 
obtained from the binned distributions of reference event samples, reweighting the Monte 
Carlo signal events with the CC03 matrix elements corresponding to various values of 
mw and T w . Two types of fits are performed for all four channels individually. In the 
first, a one-parameter fit for mw is made, where IV varies with mw according to the 
Standard Model as T w = 2.094 GeV x (m w /(80.35 GeV/c 2 )) 3 . These results provide 
the most precise value of mw- In the second, two-parameter fits are performed allowing 
mw and T w to vary as two independent parameters. Although the shape of the invariant 
mass spectra are dominated by experimental resolutions, these fits are used to test the 
validity of the SM prediction for Tw and check for any correlation between the two fitted 
parameters. Technically, the matrix element calculation assumes the Standard Model 
value for T w at a given W mass, for the coupling of electrons and their neutrinos to W 
bosons and allows the width to vary freely only in the W propagator. 

At LEP1, the Z mass was defined using a running-width scheme in the Breit-Wigner 
propagator. However, a fixed-width scheme has been employed in generating all WW 
events with K0RALW. As a result, to make both mass measurements consistent with each 
other, a positive shift of 27 MeV/c 2 is applied to the extracted W mass [29]. The 
corresponding shift to the fitted width of 0.7 MeV is not significant. 

The statistical error on mw and IV is computed from the fits to the data distributions. 
Also, a large number of subsamples are studied, each with the same number of events 
observed in the data, to evaluate the expected errors. 

The selection efficiency is found to be independent of the W mass. The variation of 
the total signal cross section with mw affects the purity of the selected events and is taken 
into account, whereas its dependence on T w is assumed to be negligible. 

The reweighting procedure was tested at 189 and 207 GeV by comparing the fitted with 
the input mass in each channel individually for four independent 4f Monte Carlo samples 
generated with m w values of 79.850, 80.100, 80.600 and 80.850 GeV/c 2 . The relationship 
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between the fitted and true masses was found to be linear for all channels over this range. 
The best straight line fits through the points are consistent with calibration curves of 
unit slope and zero bias, within the statistical precision of the test. Small deviations are 
observed in the ez/qq channel from which a systematic uncertainty is derived (Sec. 8.4). 

Table 2 gives the expected and observed numbers of events from all contributing 
processes for each channel which satisfy the kinematic fitting criteria after all window cuts 
are applied. The numbers of expected WW events are calculated with 0{a) corrections 
using the standard 4f reference samples generated at m w = 80.35 GeV/c 2 . 

Table 2: Expected numbers of events corresponding to the whole data sample (183-209 GeV) for signal 
and background processes after all selection, quality and window cuts for the four categories of events 
used in the extraction of mw and r\y All WW events are regarded as signal in the calculation of 
the quoted purities per channel. The signal cross sections are determined with mw=80.35 GcV/c 2 and 
r\v=2.094 GeV and the 0(a) correction is applied. 



Process 


4q 


ez/qq 


/iz/qq 


rz/qq 


WW — > qqqq 


4264 


0.1 


0.0 


5.0 


WW — > ez/qq 


2.1 


1217 


0.1 


120.4 


WW — »■ fiuqq 


1.9 


0.5 


1295 


41.7 


WW — > rz/qq 


10.2 


41.5 


41.6 


959.4 


qq(7) 


591 


17.9 


0.6 


35.4 


ZZ 


95 


2.2 


4.3 


23.8 


Zee 


2.2 


7.4 


0.0 


16.3 




0.0 


0.0 


0.0 


0.7 


TT 




0.2 




0.4 


77 — > TT 




0.0 




0.1 


77 — >hadrons 




0.4 




0.2 


Predicted events 


4966 


1288 


1342 


1203 


Observed events 


4861 


1259 


1371 


1226 


Purity (%) 


86.1 


97.8 


99.6 


93.6 



6.1 The qqqq channel 

The two-dimensional reweighting fits used in the previously published analyses at 183 
and 189 GeV [4, 5] are replaced by three-dimensional (3-D) fits which better exploit 
the available information from each event. The following three estimators were selected: 
(i) the 5C fitted mass, M 5C , (ii) a random choice of one of the 4C di-jet unrescaled 
masses, M 4C and (iii) the kinematic fit error on the 5C mass, o"m 5C - Using a binned 
3-D probability density function, a maximum likelihood fit is performed to the data 
within the following acceptance windows: 70 < M 5 c < 90 GeV/c 2 , < (Tm 5C < 4 
GeV/c 2 and 60 < M 4C < 110 GeV/c 2 for both the one and two-parameter fits. The 
allowed fit range for T w is loosely constrained to 1.1 < T w < 4.1 GeV. Bin sizes 
in the probability density distribution of the 5C and 4C masses are chosen for signal 
and summed backgrounds separately such that the number of events of each type per 
bin is approximately constant. The third dimension is subdivided into four bins chosen 
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dynamically to equalise the number of signal events in each bin. This binning is kept for 
the summed background. The fitted mass is extracted in each of these bins in the third 
dimension and the likelihoods combined to determine the final mass and error. To avoid 
any bias, the minimum number of signal Monte Carlo events per 3-D bin is 200. 

Fig. 3(a) shows the mass distribution from the 5C kinematic fits to the data before the 
window cuts between 70 and 90 GeV/ c 2 are applied. For comparison the mass distribution 
predicted from the simulation, reweighted to the fitted W mass in data, is superimposed. 

6.2 The ez/qq and fiuqq channels 

The following variables are used to form a three-dimensional (3-D) probability density 
function: the 2C mass, M 2 c, where the leptonic and hadronic masses are constrained 
to be equal, the kinematic fit uncertainty on the 2C mass, <Jm 2C an d the 1C hadronic 
mass, Miq. The event-by-event correlation between M^q and Mic was found to be 
43% at 189 GeV. By construction, the 3-D probability density function from the 
simulation takes into account all correlations amongst the three variables and leads to 
an improvement in statistical precision compared with a 1-D fit. Using a binned 3-D 
probability density function, a maximum likelihood fit is performed to the data within 
the following acceptance windows: 70 < M 2 c < 90 GeV/c 2 , < &m 2 c < 10 
GeV/c 2 , 60 < M^q < 110 GeV/c 2 and with the fitted Tw being constrained in the 
range 1.1 < T-w < 4.1 GeV. The bin sizes for the Monte Carlo events are chosen using 
the same criteria as for the qqqq channel. The binning of the 3-D probability density 
function has 3 intervals along the event-by-event error axis. A stable mass value and 
statistical error are obtained when the minimum number of Monte Carlo events in any 
bin is 200 or greater. 

Figs. 3(b) and (c) display the mass distributions for data resulting from the 2C 
kinematic fits to these semileptonic final states together with the predictions from the 
simulation. 

6.3 The Tisqq channel 

For ruqq candidates, a 2-D reweighting fit uses the 1-C hadronic mass, M^j and its 
uncertainty, cr M w, from the kinematic fit. The events must be within the following mass 

and error acceptance windows: 70 < M{§ < 90 GeV/c 2 and < a M w < 10 GeV/c 2 . 
In this channel, the allowed fit range for Tw is 0.9 < T\y < 4.3 GeV. The binning of 
the 2-D probability density function has four intervals along the event-by-event error axis 
and 60 intervals of varying size along the 1C mass axis. 

Figs. 3(d) displays the mass distribution resulting from the 1C kinematic fits to the 
data events together with the prediction from the simulation. 
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Figure 3: Mass distributions in the: (a) 4q, (b) ez/qq, (c) /ii^qq and (d) ruqq channels for data (points 
with error bars), non-WW background (shaded area) and signal+background Monte Carlo with mw 
values set to those fitted from each individual channel (solid line histogram). For the 4q channel, the 
distribution shows the 5C kinematically fitted dijet masses before window cuts are applied. For the ivqq 
channels, the distributions show the 2C (or 1C) kinematic fits before window cuts. 
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7 Studies on colour reconnection 



The W bosons decay at a short distance from each other (1/T « 0.1 fm), so that in the 
qqqq channel their decay products hadronise closely in space time at the typical hadronic 
scale of ~ 1 fm. An interaction between the partons from different W decays may then 
occur. 

At the perturbative level, the shift in the reconstructed W mass due to single gluon 
exchange is suppressed by the square of the number of colours and by an additional factor 
of rw/mw- The mass shift is of the order of a few MeV/ c 2 [19]. However, when the scale 
of gluon exchange is not large compared with Tw, non-perturbative colour reconnections 
(CR) in the parton cascades may lead to much larger mw shifts. 

At energies well above the pair production threshold, as in the present data set, the 
final state QED interconnection between the W's induces a shift in mw of order a em Tw /ft, 
which is a few MeV/c 2 [30] and insignificant compared with the uncertainties from non- 
perturbative QCD. 

7.1 Monte Carlo models 

At the non-perturbative level, all phenomenological implementations of CR within existing 
hadronisation models predict that the particle flow distributions per event are modified 
with the low momentum particles in the inter-jet regions being most affected. Any effect 
on high momentum particles would occur only when a jet from a W~ is aligned with 
another from a W + . Such a topology would not survive the 4-jet selection. The effect 
of CR on the fitted m w is studied using the following variants of the parton evolution 
schemes: 

(a) SKI, SKII, SKir [19] and GAL [20] in JETSET, 

(b) 2-step variants AR2 and AR20 [21] in ARIADNE and 

(c) HWCR in HERWIG [17]. 

As formulated, the SK versions in JETSET predict no effect at the Z and therefore, 
unlike the other variants, cannot be calibrated with Z data. The probability of an event 
to be reconnected depends on the string overlap between partons from the two W decays. 
In SKI, this is governed by a freely adjustable 'string' overlap parameter, k{, whereas 
the predictions of the SKII and SKII' variants are fixed once the string parameters are 
fitted in JETSET. When k\ is set to 0.65, the fraction of reconnected events is similar to 
SKII (29.2%) and SKII' (26.7%). However, SKI(A; i =0.65) predicts a larger shift in m w 
than the other SK versions. The authors state that all SK models are equally valid [31]. 

The GAL implementation within JETSET allows string rearrangements to occur by 
colour exchange with the probability for reconnection depending upon the reduction in 
total string area between the old and new configurations. After tuning at the Z on global 
event properties, the fitted value of a non-perturbative free 'strength' parameter, Rq, is 
found to be 0.04 correlated with the shower cut-off, Qo, of 1.57 GeV/c. The author 
recommends a larger value for Rq of 0.1 from fits to deep-inelastic scattering data which 
would lead to a correspondingly larger mw shift (~100 MeV/c 2 ). 

For AR2, both intra- W and inter- W reconnections are allowed between all dipoles 
with the same colour indices formed from emitted gluons with energy E g . The parton 
cascade is performed in two steps (i) allowing only intra-W reconnections with E g > Fw 
and (ii) allowing also inter-W reconnections but only for E g < r W - For AR20, no CR is 
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applied either between or within the W's. In principle, the predicted net shift in mw due 
to inter-W reconnections is determined from the difference found between AR2 and the 
corresponding variant, AR21, where only intra- W reconnections are allowed. However, in 
practice it is found that the difference between AR20 and AR21 when tuned at the Z is 
not significant. Thus, the mw shift is taken from the comparison of AR2 and AR20 fitted 
events. 

In HWCR, the criterion for allowing colour reconnections is based on the reduction in 
space-time distances within the colour singlet clusters at the end of the parton shower. 
The reconnection probability is set to 1/9 for allowed re-arrangements. The parameter 
VMIN2, the minimum squared virtuality of partons, is set to 0.1 (GeV/ c 2 ) 2 . 

Table 3 gives the predicted mass shifts <5mw = mw(CR) — mw(noCR) from these 
models averaged over CM energies from 183 to 209 GeV. Details of the parameter settings 
used in the models are given in Appendix A. 

Table 3: Predicted W mass shifts from various models averaged over all CM energies. 



Model 


5m w (MeV/c 2 ) 


SKI (Jfei=0.65) 


+39±2 


SKI (jfei=1.0) 


+56±2 


SKII 


+6±8 


SKII' 


+4±8 


AR2-AR20 


+54±5 


HWCR 


+39±4 


GAL (i?o=0.04) 


+44±8 



The predicted mass shifts from the models tunable at the Z: GAL, AR2 and HWCR, 
range from 40 to 55 MeV/c 2 , suggesting for consistency that the value of k\ in the SKI 
model should be of order 0.8. 

To examine the validity of some of the tunable models of CR at the Z, the particle 
distributions in selected three-jet events were compared specifically with the predictions 
of AR2 and GAL [32] . If it can be assumed that the behaviour of colour rearrangements 
in the parton cascades of Z decays is the same as for WW, these observations suggest 
that the two models overestimate the effects on mw from CR. 

7.2 Data Analysis 

Keeping the originally reconstructed jets in each selected event, the W mass analysis is 
repeated twice, either removing all low momentum particles (PCUT analysis) [27], or 
rejecting particles outside cones directed along the four jet axes (CONE analysis) [33]. 
The difference from the mass measured without these additional cuts, called the standard 
analysis, is a sensitive observable of the CR effect according to all the above models. 

For each of five values of the particle momentum cut off from 1 to 3 GeV/c in the 
PCUT analysis, each jet energy and angle is recomputed. In the CONE analysis, each jet 
energy is kept unchanged, whilst its three-momentum is recomputed from the vector sum 
of its remaining participating particles, rescaled by the ratio of the original jet energy 
to the energy of the particles inside the cone. Seven values of the cone opening angle R 
are used from 0.4 to 0.9 radians. Studies show that fragmentation uncertainties increase 
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rapidly for momentum cut-offs beyond 3 GeV/ c or cone angles smaller than 0.4 radians. 
These values were found to provide optimal balance between statistical and systematic 
uncertainties on my/. 

Figure 4 shows the expected variation of the mass due to CR as a function of the cut 
for the tuned AR2, HWCR and GAL models in the 183 to 209 GeV energy range. The SKI 
predictions for two values of k\ are also included. The predictions for each of the eight CM 
energies are combined using the relative integrated luminosities of the data. Table 4 lists 
the m w shifts, 5m w PCUT and 5m w CONE , for the PCUT(=3 GeV/c) and CONE (R=0.4 
rad) reconstructions respectively. The corresponding mw shifts in the standard analysis 
5mw° are shown for comparison. Within errors the shifts for each reconstruction are 
comparable for all tuned models and consistent with SKI (ki—1). 



Table 4: Predicted W mass shifts ((5mw) from various models averaged over all CM energies for the 
CONE (R=0.4 rad) and PCUT (=3 GeV/c) reconstructions (units in MeV/c 2 ). 



Model 


5m w ° 






SKI (M=l-0) 


+56±2 


+19±4 


+23±3 


AR2-AR20 
HWCR 

GAL (#0=0.04) 


+54±5 
+39±4 
+44±8 


+ 17±8 
+13±7 
+27±12 


+20±6 
+14±7 
+22±11 
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Figure 4: 5m w versus (a) PCUT in GeV/c and (b) inverse CONE radius (R) in rad" 1 for SKI (2 h 
values), AR2, HWCR and GAL models in the qqqq channel. 



For the data collected at all CM energies combined, Fig. 5 shows the mass difference 
Am-yv between a PCUT or CONE reconstruction and the standard mass analysis. The 
slopes are fitted with the full correlation matrix included and amount to —11 ± 16 
(MeV/c 2 )/(GeV/c) for the PCUT analysis and +9±19 (MeV/c 2 )/^- 1 ) for the CONE. 
Both values are compatible with no effect. 
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Figure 5: Am w versus (a) PCUT in GeV/c and (b) inverse CONE radius (R) in rad -1 for qqqq data. 
The dashed lines denote the fitted slopes and ±lcr errors. The correlation with respect to the standard 
analysis is taken into account in the error on the mass difference for each reconstruction. 



A cross check was performed on all the semileptonic channels where no CR effect 
between the decay products of the different W's can be present. The mass analyses in 
the ez/qq and /iz/qq channels were repeated for PCUT and CONE following the same 
kinematic fit procedure as used in the rz/qq channel where only the hadronic jets are 
included. Figure 6 shows the corresponding mass differences for each cut value relative 
to the standard analysis after combining the results statistically from the ez/qq, /iz/qq 
and rz/qq channels. No significant instability is observed. The combined iuqq channels 
represent a sample of size similar to the size of the qqqq channel and give a slope of 
+4 ± 21 (MeV/c 2 )/(GeV/c) for the PCUT analysis and +16 ± 27 (MeV/c 2 )/^- 1 ) for 
the CONE, which are not significantly different from zero. 

A limit on 5m\y can be inferred from a comparison between the slopes observed in the 
data and those from the CR models. For each model, pseudo-data samples were built, 
combining all the CM energy points weighted by their respective integrated luminosities. 
In the case of the SKI model, 20 different values of the k\ parameter are chosen, ranging 
from to 100. 

The SKI model predicts a clear correlation between the mass shift for the standard 
reconstruction, <5m\y°, and the slope of the mass difference as a function of the PCUT or 
CONE cuts as shown in Fig. 7. The clustering of the slope values from AR2, HWCR and 
GAL, around -10(MeV/c 2 )/(GeV/c) for PCUT and similarly for CONE corresponds to 
the previously described values of <5mw quantified in Table 4. 

The covariance between the slopes, S cone and S pcvLt , is computed as well as the 
resolution on the slopes from a Gaussian fit using the pseudo-data samples at each ki. The 
average correlation between the PCUT and CONE slopes is 51 % with little dependence 
on k Y . The RMS errors on the slopes are 18 (MeV/c 2 )/(GeV/c) for the PCUT analysis 
and 19 (MeV/ c 2 )/( rad -1 ) for the CONE analysis, in agreement with the values obtained 
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Figure 6: Am w versus (a) PCUT in GeV/c and (b) inverse CONE radius (R) in rad 1 for data from 
the e,/i,r^qq channels combined, fitting with the jets alone in each case 
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Figure 7: Slope of the mass difference relative to the standard analysis as a function of Smw° f° r 
(a) the PCUT and (b) the CONE reconstructions. The dashed line is a straight line fit to the SKI 
points represented by white circles. The black symbols represent AR2 (circle), HWCR (square) and GAL 
(triangle) predictions. The horizontal bands represent the measured slopes with their ±1ct errors. 



in the data sample. A x 2 fit, defined as follows: 

£(^r a - S^ c {x))a-j{Sp a - Sf°{x)) 

a/3 

20 



where a,f3 signify PCUT and CONE respectively and a a /3 the covariance matrix, is used 
to extract the 68% CL Gaussian upper limit on x. The parameter x can be either k\ or 

The x 2 curve is shown in Fig. 8 as a function of 5m^°. The upper limit on Sm^° has 
been set to +78 MeV/ c 2 , corresponding to the value at which the integral of the Gaussian 
likelihood from zero is 68% of the full integral over the allowed (positive) range. 
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Figure 8: A\ 2 as a function of the mass shift Smyj in the standard analysis. The dotted line is 
with statistical errors only, the dashed line includes fragmentation errors, and the full line includes all 
systematic uncertainties. 



The fragmentation uncertainty on this limit is estimated from the average bias on 
the slopes between an ARIADNE or HERWIG sample and the JETSET reference sample, 
repeated over the pseudo-data samples. The slope biases from ARIADNE, being largest, are 
used as diagonal terms of a systematic error matrix with 100% correlation between the 
PCUT and CONE slopes. This matrix is added to the statistical error matrix resulting 
in an increase in the mass limit to +86 MeV/ c 2 . An estimate of the uncertainty in 
evaluating the statistical error matrix yields a further small increase in the mass limit 
to +87 MeV/ c 2 (corresponding to ki (68% U.L.) = 1.88). Figure 8 shows the progressive 
effect of adding these systematic uncertainties to the \ 2 fits. 

The same procedure is used to determine the mass shifts 5mw PCUT and 5m\y E - 
The upper limits on these mass shifts are given in Table 5 together with the limit from 
the standard analysis. 

Both AR2 and HWCR agree well with the SKI prediction of the slopes as a function 
of Smyj. However, the actual limit from SKI is larger than the AR2, HWCR and GAL 
absolute predictions and therefore is taken as the systematic uncertainty due to colour 
reconnection for each reconstruction. 
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In practice, these limits on the mass bias depend linearly on the CM energy within 
SKI. From 183 to 209 GeV, the limit varies from 45 to 105 MeV/ c 2 when no cut is applied, 
12 to 36 MeV/c 2 for the PCUT and 18 to 45 MeV/c 2 for the CONE respectively. 

Table 5: CR mass shifts (68 % C.L. upper limits in MeV/c 2 ) for the three reconstructions: standard, 
PCUT and CONE using SKI derived from (a) the purely statistical analysis and (b) incorporating all 
systematic uncertainties. 





(a)SKI(stat.) 


(b)SKI(stat.+syst.) 


5m w ° 


78 


87 


5m w PCUT 


25 


27 


Sm w CONE 


32 


35 



8 Systematic uncertainties 

Systematic uncertainties in the measurement of to-w and Tw arise from an incomplete 
description of the WW production process, inaccuracies in the simulation of event 
reconstruction in the detector and the modelling of the W decays to di-jets. The following 
subsections describe all the systematic uncertainties evaluated for the standard analysis 
in each of the four event categories. They are also determined in the qqqq channel for the 
CONE (R=0.4 rad) and PCUT (=3 GeV/ c) reconstructions where the potential effects 
of colour reconnection (CR) are minimised. 

The LEP energy uncertainties with year-to-year correlations are taken from Ref. [34]. 
All other uncertainties in the analysis are evaluated at 189 and 207 GeV, simultaneously 
in mw and Tw from the two-parameter fits. When combining all the measurements, any 
variation over this energy range is taken into account using a linear interpolation for the 
intermediate CM energies. Table 8 lists all the systematic uncertainties for the standard 
analysis as well as the optimal PCUT and CONE reconstructions in the qqqq channel. 
The CR uncertainty in this channel is taken into account at each CM energy. Table 9 
lists all the systematic uncertainties in the standard analyses of the three semileptonic 
channels. 

8.1 Detector simulation 

The systematic uncertainties in the detector simulation for the qqqq events are those 
arising from the quantitative comparison of the reconstructed jet four-momenta with 
the data as described in Sec. 3. For the ez/qq and /jvqq channels, the uncertainties in 
the lepton four-momenta are included and combined in quadrature with those from the 
jets. Subsidiary studies of particles within the jets have been made by comparing the 
simulation with data for the effect of photon energy miscalibration and charged hadron 
tracking discrepancies. These uncertainties are already taken into account in those quoted 
for the jets. Each uncertainty is evaluated by first comparing the mean fitted parameters 
from special pseudo-data samples with corresponding normal samples each of the size of 
the data. The mean shifts found in mw and Tw are then rescaled to correspond to the 
residual discrepancies found between data and simulation after any corrections have been 
applied (Sec. 3.3). 
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8.1.1 Isolated lepton reconstruction in ez/qq and /xz/qq events 

Specific studies (Sec. 3.3.2) have been performed for electrons and muons. In the ez/qq 
channel, the uncertainty is determined from the error (0.04%) in applying the global 
momentum correction of 0.45% combined with the percentage error of 0.0008% per GeV 
in the evaluation of the momentum scale linearity, taking into account the correlation 
(+0.78). The small biases found as a function of polar angle have a negligible effect. 

For the /zz/qq channel, the momentum uncertainty is derived from the full effect of the 
uncorrected global offset of 0.08%. The percentage error of 0.0025% per GeV in the slope 
is added in quadrature, taking into account the correlation (-0.22). 

Averaged over polar angle, the lepton momentum resolutions in the simulation are 
degraded by 13% and 8% for the electrons and muons respectively to match the data 
when averaged over all momenta. For mf, the effect is relatively insignificant but on Tw 
it is the dominating contribution to the total systematic uncertainty in each channel. 

A possible bias in the measurement of the lepton direction in the ez/qq and /iz/qq 
channels was studied by comparing the lepton track 9 and angles as measured by the 
VDET and the ITC + TPC separately [5]. No difference greater than a fraction of a mrad 
was observed. Owing to small offsets in the drift time of the TPC, the z-component of 
momentum can be biased for tracks away from 90° to the beam axis. Conservatively, the 
effect on the lepton polar angle is parametrised maximally as 2.0 x sin2#i cp ton mrad with 
respect to the beam axis. Events are generated accordingly, whilst keeping the lepton 
energy and the total momentum of the event conserved. The shift is negligible for both 
raw and r^. Any effect from possible lepton <fi angle biases is also negligible. 

Comparing the VDET and ITC + TPC track measurements [5], the spread of the 
differences in polar angle measurement for the electrons and muons was found to be of 
order 0.5 mrad. No mean discrepancy greater than 0.3 mrad between the data and Monte 
Carlo distributions was observed. Conservatively, an additional 0.5 mrad smearing has 
been applied to the simulation to compute the uncertainties attributable to the simulation 
of angular resolution. The shifts in m w are found to be negligible. 

8.1.2 Jet energy corrections before the kinematic fit 

As described in Sec. 3.3.4, the simulation of jet energies from di-jet events produced at 
the Z was compared with data in the range 30 to 70 GeV. It was shown that (a) the bias 
in the relative global energy scale does not exceed 0.5% in the central region rising to 
2.5% at low angles and (b) the relative slope of the data to simulation in the jet energy 
scale as a function of E- Jet is flat setting a limit of ±0.8 x 10~ 4 per GeV. Studies with 
special simulated samples show that the global bias has no significant effect on m\y and 
Tw for all channels. The systematic uncertainties assigned from the limit on the slope 
variation with E- )et is combined in quadrature with the shifts obtained from a mean 1% 
global discrepancy between barrel and endcaps. Disregarding the presence of b jets in the 
Z samples, introduces a shift of 0.25% in the relative global energy scale. Since this is 
only half of the bias taken into account, the effect due to b jets is negligible. 

8.1.3 Jet energy resolution 

As described in Sec. 3.3.4, the data and Monte Carlo resolutions in each cos#j Ct bin as 
determined from the RMS spread of jet energies agree to within ±2% for di-jet events at 



23 



the Z. Special samples were made where the jet four-momenta and energies are smeared 
degrading the resolution by 10% with respect to the nominal values computed from the 
kinematic fit parametrisations. The shifts found are rescaled to correspond with the 
measured difference. 

8.1.4 Jet angular bias 

Possible discrepancies in the determination of 9- iCt were studied [5] by comparing, between 
data and simulation, the direction of the charged and neutral jet components in Z di- 
jet events. The tracking detectors and the ECAL were aligned independently but high 
statistics studies performed at 91.2 GeV show that their relative polar angle alignment is 
about 1 mrad. In order to measure any angular distortions, separately constructed charged 
and neutral components of jets are selected and their polar angle directions compared in 
bins of 5 degrees in (Charged + #neutrais)/2. The simulation of the jet components is in 
good agreement with the data except in the overlap region between the barrel and endcap 
calorimeters where the difference is up to 2 mrad. The difference A(# cha r g cd — ^neutrals), 
parametrised as 0.7(0.2) - 2.4(0.6)cos# sin# mrad (errors in brackets), gives the best fit to 
these discrepancies. Further studies confirm that the effect of the global offset of 0.7 mrad 
is negligible since the uncertainties effectively cancel in cos 8j et by symmetry. The resulting 
systematic uncertainties are evaluated by applying this parametrisation without the offset 
to special Monte Carlo WW event samples assuming that A#j Ct = (Charged _ ^neutrals) /2. 

8.1.5 Jet angular resolution 

Selected di-jet events from the Z calibration runs have been used to measure the jet 
angular resolution for 45 GeV jets from the distribution of the opening angles between 
the two jets. The PCUT and CONE criteria are also applied to the jets to measure the 
variation in jet angle resolution for these reconstructions. Special event samples with 
modified resolutions, which match the data, are used to estimate the effect on m-w and 
T-w for all channels. The resulting uncertainties in m-w are very small. 

8.1.6 Jet boosts 

The accuracy of the Monte Carlo reconstructed jet masses in each channel depends 
sensitively on the simulation of the charged and neutral particle momenta and multiplicity 
distributions within the jets. Jet boosts, /Jj e t7jet, are chosen rather than masses to 
compare data with simulation since any momentum discrepancies are factored out and 
double counting minimised. Figure 9 compares the data and Monte Carlo distributions of 
log(/9j et 7jet) for jets built as in the standard analysis, integrated over all polar angles from 
(a) high statistics hadronic Z decays where the average jet momenta are close to those in 
W decays, (b) higher energy di-jets, (c) hadronic Z decays for PCUT and (d) hadronic Z 
decays for CONE. These jet samples are studied rather than those from the selected W 
pairs to avoid the possible influence of final state interactions and to benefit from high 
statistics. The study includes b-depleted samples and jets from radiative returns to the 
Z peak (Sec. 9). Table 6 gives the largest shifts obtained between data and simulation 
expressed as A log(/9j e t7jet)- The small differences between central and forward regions of 
the detector are not statistically significant. The systematic uncertainties in mw and Tw 
are derived using special Monte Carlo samples rescaled to match the biases in this table. 
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Figure 9: Distributions of jet boosts (log/3j e t7j e t) for data (circles) and MC (histogram): (a) from 
Z — > qq events (1998-2000) and (b) from high energy di-jet events (183-209 GeV) using the Durham 
jet reconstruction in the standard analysis, (c) and (d) from Z — > qq events using PCUT and CONE 
reconstructions respectively. 



8.2 Fragmentation of the W — > qq decays to hadrons 

In the previous analysis at 189 GeV [5], the uncertainty due to the modelling was 
determined mainly from the comparison of mw and r\y values using event samples in 
which fragmentation is simulated with HERWIG [17] or ARIADNE [18] in place of JETSET. A 
large uncertainty of ~35 MeV/c 2 , fully correlated between channels was assigned. It has 
been found that the variation in baryon content between the models is largely responsible. 
The baryon multiplicities predicted by JETSET and ARIADNE agree with data at the Z [35] 
whereas HERWIG generates (~ 0.5) fewer baryons per event. 
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Table 6: Largest (data-MC) shifts, Alog(/?j e t7j e t) m percent. The shifts are tabulated for the central 
region of the detector (| cos0j o t| < 0.7), the forward region (| cos#j ot | > 0.7) and both combined (errors 
are shown in brackets). 



Reconstruction 


central 


forward 


combined 


standard 

PCUT 

CONE 


0.9 (0.2) 
2.7 (0.2) 
2.1 (0.4) 


0.7 (0.2) 
1.9 (0.3) 
2.1 (0.5) 


0.8 (0.1) 
2.4 (0.2) 
2.1 (0.3) 



The uncertainties in m w for each channel in the standard analysis are reassessed 
after correcting for this effect. In the qqqq channel, the bias in m w is found to depend 
linearly on the number of protons and neutrons per event. Taking samples with 0, 2, 
4, 6 and 8 nucleons per event, the slope of the bias for all three models is statistically 
equivalent and found to be 20.1±0.8 MeV/c 2 per nucleon. A similar linear behaviour is 
seen in the ez/qq, /xz/qq, and rz/qq channels. The W mass differences between the models 
due to the variation in their baryon content is evaluated from the linear dependences 
in each channel assuming that they apply over the entire range of baryon multiplicities. 
For HERWIG-JETSET and ARIADNE- JETSET , the mass shifts before and after correcting for 
the differences in baryon content are given in Table 7. After correction, the differences 

Table 7: For the standard analysis in the 183-209 GeV range, the mean W mass differences between 
MC samples of HERWIG and ARIADNE relative to JETSET are tabulated for each channel before and after 
correcting for the difference in baryon content. 





HW-JT 


AR-JT 




Uncorrected 


Corrected 


Uncorrected 


Corrected 


qqqq 


+ 12±8 


-7±8 


+3±9 


+5±9 


ez/qq 


+25±8 


+3±8 


+ 1±9 


+6±9 




+10±8 


-8±7 


-11±8 


-7±8 


rz/qq 


+40±11 


+ 15±11 


+5±13 


+6±12 



between HERWIG and JETSET become insignificant. All three fragmentation models agree 
within statistical error for all channels. The systematic uncertainty is set to 10 MeV/ c 2 for 
the standard analysis, coherent in all channels. The variation in baryon content between 
the models has no significant effect on the values fitted for Tw- 

For the PCUT and CONE reconstructions in the qqqq channel, the uncertainties are 
determined from comparing the same event samples simulated with ARIADNE and JETSET 
where the variation in baryon content is not significant. Any residual differences (AR- 
JT) are due to other effects unrelated to baryon multiplicities and are taken to represent 
the uncertainties for these reconstructions. The differences are larger than the standard 
analysis but comparable for both reconstructions. 

8.3 Radiative corrections 

The uncertainties in the theoretical treatment of QED initial state radiation (ISR) and 
Coulomb corrections in K0RALW as well as Next-to-Leading 0(a) corrections in YFSWW3 
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are determined for each channel and reconstruction by comparing Monte Carlo samples 
with appropriate event weighting. The estimated shifts from each study are combined in 
quadrature. 

8.3.1 Missing ISR corrections 

Initial state radiation is simulated in KORALW up to 0(a 3 L 3 ), i.e. up to third order in the 
leading-log approximation. The effect of missing higher order ISR terms beyond 0(a 3 L 3 ) 
on the measurement of m w and T w is estimated by measuring the respective shifts when 
this QED computation is downgraded to 0(a 2 L 2 ) as originally suggested in Ref. [36]. Each 
event in a specially generated KORALW sample is weighted according to the calculated ratio 
of second to third order squared matrix elements: 0(a 2 L 2 )/0(a 3 L 3 ). Treated as data, 
fits are made to the weighted events selected in each channel and compared with those 
from the corresponding unweighted events to evaluate the shifts. The shifts in mw and 
T w are less than 1 MeV in all channels. 

8.3.2 Coulomb corrections 

The unweighted events from KORALW include non-factorizable QED corrections [37] 
which effectively "screen" the Coulomb interaction [30] between the two W's. It is 
suggested [38] that the difference between this "screened" Coulomb correction and no 
Coulomb correction can be used to assess an uncertainty. The differences are found to be 
less than 3 MeV in all channels for both m w and T w . 

8.3.3 Next-to-Leading 0(a) corrections 

These corrections are large, ranging for mw from ~10 MeV/ c 2 in the qqqq channel to ~20 
MeV/ c 2 in the ez/qq channel. Studies have shown [39] that the theoretical implementation 
of these corrections in RacoonWW [40] are in good agreement with YFSWW3. The following 
two possible contributions to the uncertainties in these corrections using YFSWW3 are 
considered. 

(i) The main effect of the NL 0(a) corrections is to modify the W final state radiation 
(WSR) pattern of photons. In YFSWW3, the infra-red (IR) contributions to WSR and WSR- 
ISR interference are exponentiated to infinite order including non-IR Next-to-Leading 
contributions. The uncertainty in this calculation is estimated as suggested in Ref. [38] 
by removing the additional non-IR contributions. The effect is found to be small, less 
than 2 MeV/c 2 in m w for all channels and reconstructions. The shifts in T w are similar. 

(ii) In calculating the weight per event from the YFSWW3 program, the recommended 
recipe by the authors [9] is an additive correction in which the Double-Pole approximation 
(DPA) for doubly resonant W's is applied only to the CC03 part of the event weight. 
An alternative recipe would be to apply the NL correction also to the difference 
between the 4f and CC03 contributions - the so-called multiplicative NL correction. The 
additive correction is chosen as the default. The systematic uncertainties in mw and 
Tw are estimated by measuring the difference between the additive and multiplicative 
implementations. These differences are less than 1.5 MeV/c 2 in all channels. 
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8.4 Calibration curves 

As stated in Sec. 6, the reweighting procedure was tested by comparing the fitted with 
the input W masses and widths in each channel individually. No deviations were observed 
in the fitted slopes or intercepts of the produced calibration curves except in the euqq 
channel. Combining statistically the fitted masses at five points between 79.85 and 80.85 
GeV/ c 2 from 189 and 207 GeV pseudo-data in this channel, the calibration curve for mw 
is found to be linear but with a slope of 0.954±0.023. At the measured mass, this deviation 
from unity corresponds to an uncertainty of 10 MeV/ c 2 in mw- A similar analysis for T\y 
found no significant effect in the ez/qq channel. An upper limit of 20 MeV is assigned as 
the systematic uncertainty from the statistical precision of the test. 

8.5 Background contamination 

The expected numbers of events in each channel included in the reweighting fits from 
non-WW background processes are shown in Table 2. 

The dominant background in the qqqq channel is qq(7) (14% of all selected 
events) followed by ZZ (2%). The normalisations of the these contributions are varied 
conservatively by 5% and 10% respectively and the consequent shifts added in quadrature. 
In addition, the uncertainty in the fragmentation modelling of the qq(7) events is 
estimated by replacing JETSET with ARIADNE; its impact is significant only for m w . The 
effect of any qq hadronisation uncertainty in the ZZ contribution is very small and has 
been ignored. 

In all ivqq channels the contamination is relatively small but also dominated by events 
from the qq(7) and ZZ processes. Their rates are also varied by 5% and 10% respectively 
to produce the quoted uncertainties in mw and T W - Any effect from hadronisation is 
found to be insignificant. 

For all channels, the Zee contributions are flat in the defined mass windows and their 
effects on m w and T w are negligible. 

8.6 Final State Interactions in the 4q channel 
8.6.1 Colour reconnection 

The studies on the mass shift coming from possible colour reconnection between decay 
products of the W pairs have been discussed in section 7. 

The mass differences obtained when comparing PCUT or CONE analyses with the 
standard analysis give no indication of an effect within our data statistics; nor do the 
differences in IV (Fig. 10). The upper limit derived with SKI exceeds the predictions 
of the HWCR, AR2 and GAL models. The predictions of AR2 and GAL on particle 
distributions in three jet events at the Z have been studied in Ref. [32]. They are 
disfavoured by the data. Disregarding the SKII and SKIL models which predict a very 
small effect, the 68% upper limits on <5mw and STyj obtained with the SKI model are 
taken as conservative estimates for the systematic uncertainties from CR. The averaged 
values of 5m-w in Table 5 are determined using only the statistical errors in m w at each 
CM energy. In the final optimising combination procedure, the systematic uncertainties 
for CR quoted in Table 8 are computed taking into account both the statistical and 
systematic uncertainties from all sources. For the standard analysis where the CR 
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uncertainty dominates and varies significantly with CM energy, the value quoted in Table 
8 is consequently reduced from 87 to 79 MeV/ c 2 . 
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Figure 10: W width differences A T w versus (a) PCUT and (b) inverse CONE radius for data from 
the qqqq channel. The errors take into account the correlation at each point with the standard analysis. 



8.6.2 Bose Einstein correlations 

The presence of Bose-Einstein correlations between the decay products of the two W's 
in the WW — > qqqq selected events could influence the W mass measurement [22, 36]. 
When simulated events are modified according to the JETSET-LUB0EI model [23] of Bose- 
Einstein correlations between the Ws, tuned on hadronic Z decay data, a shift on mf of 
— 32±5 MeV/ c 2 is predicted in the standard analysis. This shift is reduced in the optimal 
CONE or PCUT analysis by a factor of two. The ALEPH dedicated analysis of Bose- 
Einstein correlations based on the comparison of like-sign and unlike-sign pion pairs using 
the so-called "mixed" method, is described in Refs. [41, 42]. The data are in agreement 
with the hypothesis where Bose-Einstein correlations are present only for pions coming 
from the same W. The JETSET-LUB0EI model with Bose-Einstein correlations applied 
also on pions from different W bosons is disfavoured by up to 4.7a using the different 
variables studied. The systematic uncertainty on is determined from the fraction 
of the full prediction of this model which is consistent with these experimental results, 
using the value predicted with and without Bose-Einstein correlations between pions from 
different Ws. This fraction is — 5%±22%, giving an uncertainty on of 6 MeV/c 2 , if 
a linear dependence between the mw shift and the value of this fraction is assumed. 



8.7 LEP energy 

The LEP beam energies were recorded every 15 minutes, or more frequently if required by 
the machine conditions. The instantaneous values recorded nearest in time to the selected 
events are used in the analysis. For the year 2000, as the CM energy was continuously 
increased, the dataset is split into two samples, the first integrating data at energies from 



29 



202.5 GeV to 205.5 GeV centred at 204.86 GeV and the second including all data above 
205.5 GeV centred at 206.53 GeV. The effect on m w of any discrepancy between the 
data and generated reference beam energies was investigated and found to range from 8 
MeV/c 2 per GeV difference at 189 GeV to 16 MeV/c 2 per GeV at 207 GeV. The resulting 
uncertainties at each CM energy are small compared with the LEP energy uncertainties 
and have been ignored. 

The year-on-year correlated uncertainties in the LEP beam energy E LEP taken from 
Ref. [34], are used to determine the quoted systematic uncertainties in mw and r w . The 
relative uncertainty on mw for the rz/qq and qqqq channels is obtained directly from the 
relative error in -Elep whereas for the ez/qq and /iz/qq channels, the relative uncertainty 
is 0.9x A_Elep/-Elep- The effect of smearing in the event-by-event collision energy [34], 
which also introduces a longitudinal boost in the CM frame of ALEPH, are both taken 
into account in the evaluation of the uncertainty in r W - 



Table 8: Summary of the systematic errors on tow and % averaged over 183-209 GeV in the qqqq 
channel for the standard, PCUT (= 3.0 GeV/c) and CONE (R=0.4) reconstructions. 





Am w (MeV/c 2 ) 


Ar w (MeV) 


Source 


standard 


PCUT 


CONE 


standard 


PCUT 


CONE 


Jet energy scale/linearity 


2 


2 


3 


2 


12 


4 


Jet energy resoln 





1 





7 


9 


10 


Jet angle 


6 


6 


6 


1 


3 


3 


Jet angle resoln 


1 


3 


2 


15 


18 


9 


Jet boost 


14 


15 


11 


5 


5 


4 


Fragmentation 


10 


20 


20 


20 


40 


40 


Radiative Corrections 


2 


2 


2 


5 


7 


7 


LEP energy 


9 


10 


10 


7 


7 


7 


Ref MC Statistics 


2 


3 


3 


5 


7 


7 


Bkgnd contamination 


8 


5 


5 


29 


31 


32 


Colour reconnection 


79 


28 


36 


104 


24 


45 


Bose-Einstein effects 


6 


2 


3 


20 


10 


10 
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Table 9: Summary of the systematic errors on raw and T\y in the standard analysis averaged over 
183-209 GeV for all semileptonic channels. The column labelled iisqq lists the uncertainties in m\y used 
in combining the semileptonic channels. 





Am w (MeV/c 2 ) 


Ar w (MeV) 


— 

Source 


ez/qq 


/iz/qq 


rz/qq 


Ivqq 


ez/qq 


/iz/qq 


rz/qq 


£vqq 


e+fj, momentum 


3 


8 


- 


4 


5 


4 


- 


4 


e+/i momentum resoln 


7 


4 


- 


4 


65 


55 


- 


50 


Jet energy scale/ linearity 


5 


5 


9 


6 


4 


4 


16 


6 


Jet energy resoln 


4 


2 


8 


4 


20 


18 


36 


22 


Jet angle 


5 


5 


4 


5 


2 


2 


3 


2 


Jet angle resoln 


3 


2 


3 


3 


6 


7 


8 


7 


Jet boost 


17 


17 


20 


17 


3 


3 


3 


3 


Fragmentation 


10 


10 


15 


11 


22 


23 


37 


25 


Radiative corrections 


3 


2 


3 


3 


3 


2 


2 


2 


LEP energy 


9 


9 


10 


9 


7 


7 


10 


8 


Calibration (ez/qq only) 


10 






4 


20 






9 


Ref MC Statistics 


3 


3 


5 


2 


7 


7 


10 


5 


Bkgnd contamination 


3 


1 


6 


2 


5 


4 


19 


7 



9 Radiative returns to the Z peak 

Radiative events e + e~ — > ffy where the invariant mass of the ff system is in the vicinity of 
the Z mass are selected over the full CM energy range \fs = 183-209 GeV. The hadronic 
final states producing two jets are analysed using the same jet reconstruction methods as 
applied to the W — > qq decays, providing a cross check of the W mass reconstruction. 
Furthermore, the analysis of the /i + /i~7 channel provides a direct measurement of the 
LEP energy [43] reaching an interesting precision when combined with the other LEP 
experiments [44]. 

Candidate qqy events are required to have at least eight good tracks with total energy 
exceeding 20% of the nominal CM energy. The scalar sum of the transverse momentum 
components of the good tracks must further exceed 12% of the nominal CM energy. 
Identified photons with energy exceeding 5% of the nominal CM energy, and isolated from 
the good tracks, are rejected and ignored in the analysis. As described in Section 3.2 for 
the W mass analysis, all energy flow objects below 15° to the beams are rejected, and the 
same thresholds applied to ECAL and HCAL neutral objects. Reconstruction efficiencies 
for good tracks in data have been compared with the simulation at y / i=91.2 GeV revealing 
lower efficiencies for soft tracks from data in the forward direction. Correction factors have 
been applied to the simulation for tracks with | cos#| > 0.6 and pp < 5 GeV. Both this 
last correction for forward good tracks and the removal of neutral objects near the beam 
line are of crucial importance for the correct simulation of the forward region and the 
following reconstruction of the hadronic Z mass. The forward tracks correction was not 
applied to the generated WW events. However, its effect was found to be closely correlated 
to the jet angular bias discussed in Sec. 8.1.4 and covered by the corresponding systematic 
uncertainty. Figure 11 shows the distributions of the jet boosts in log(/3j e t7jet) comparing 
data and simulation separately for central and forward reconstructed jets. 
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As in previous ALEPH studies of hadronic radiative returns [43], and similarly to the 
W mass reconstruction, the Z mass is obtained by clustering the hadronic system into two 
jets with the DURHAM-PE algorithm, and performing a kinematic reconstruction based on 
fixing the jet velocities to their measured values but rescaling their energies to conserve 
four-momentum. It is assumed that the ISR photon is emitted along the beam line, and 
thus the boost of the produced qq system is in the opposite direction. The di-jet rescaled 
Z mass can be expressed in terms of the polar angles (9 1: 9 2 ) and velocities (/3i, /3 2 ) of the 
two jets as 



M 2 z = s 



fii sin 9\ 


+ ft 


sin 9 2 - P1P2 


sin(0i + 2 )| 


f3\ sin 6*i 


+ P2 


sin 9 2 + PiP 2 \ 


sin(0i + 2 )| 



Requiring a di-jet rescaled mass in the window 75< M z < 115 GeV/c 2 , a total of 25908 
events are selected from the data compared with 25904 events from the simulation. The 
expected signal purity is 93.8%. 

The shift of the Z mass peak is measured by means of an unbinned likelihood fit 
to a p.d.f. built from reference distributions. The calibration of the fit is done with 
pseudo-data samples and the small bias is corrected. 

Various sources of systematic uncertainties on the AM Z measurement have been 
considered. Background uncertainties have been evaluated by varying the expected 
contribution of the background events as Zee (±50%), Wez/ (±25%) and 77qq (±100%), 
leading to a combined effect on AM Z of 16 MeV/ c 2 . Fragmentation systematics have 
been evaluated by comparing results obtained with different models JETSET, ARIADNE 
and HERWIG leading to an uncertainty on AM Z of 19 MeV/c 2 . For the calorimeter 
systematics that have an impact on the jet boost, different shower simulations have 
been used (Section 3.1) and in particular the use of FULLSIM leads to a difference in 
AM Z of 30 MeV/c 2 . For the tracking affecting the jet angles, half of the full effect 
due to reconstruction inefficiencies in the forward direction is taken as the systematic 
uncertainty of 16 MeV/ c 2 . The uncertainty related to the ISR model is estimated to be 7 
MeV/ c 2 . The uncertainty coming from limited Monte Carlo statistics is dominated by the 
fit calibration uncertainty and is 12 MeV/c 2 . Using different fit methods the uncertainty 
due to the fit method is estimated to be 20 MeV/c 2 . Possible global biases of 0.2 mrad 
on the polar angle measurements of the good tracks lead to an uncertainty of 24 MeV/ c 2 . 
The combined systematic uncertainty on AM Z due to all the above sources is then 54 
MeV/c 2 . 

The resulting shift in the di-jet Z mass peak reconstruction in radiative events is 

Am z = ±40 ± 30(stat.) ± 54(syst.) MeV/c 2 , 

which is consistent with zero. This conclusion remains unchanged when the di-jet rescaled 
Z mass is evaluated using jets built with CONE, PCUT or with good tracks only. The 
jet reconstruction methods studied here are applied to the determination of the W mass 
and similar uncertainties are used. Thus, these conclusions give further confidence in the 
W mass analysis. 

The previous result was obtained using the beam energies supplied by LEP. If the 
measurement is in turn interpreted as a shift of the nominal LEP2 CM energy, where m z 
is kept fixed to the published value, 

Av^ = -86 ± 64(stat.) ± 116(syst.) MeV 
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Figure 11: Distributions of jet boosts (log/3j e t7j e t) for hadronic events at 183-209 GeV comparing Data 
and simulation (a) in the central region (cos(#j e t) < 0.7) and (b) in the forward region (0.7 < cos(#j Ct ) < 
1.0). 



is obtained, which again is in good agreement with zero. 

The analysis of muon pairs from the process e + e~ — > Z7 — > fi + fi~ / y provides an 
additional check on the reconstruction of the LEP beam energies. Two variables are used 
(i) the plain invariant mass Mi 2 , defined as M 2 2 = 2PiP 2 (l — cos #12), where P\,Pi are 
the momenta of the two muons, 812 is the angle between them and (ii) the angular mass 
m\i given by 

2 sin 81 + sin 82 — I sin 8 12 1 
12 sin 81 + sin 82 + | sin 8i 2 \ 

Selected di-muon events are required to be in the range M12 > 60 GeV/c 2 and 
80 < 77112 < 100 GeV/c 2 . A total of 976 events are selected from the data and 971.2 
are expected from the simulation, with an expected signal purity of 93.4%. 

Any discrepancy between data and simulation in the M 12 and m 12 distributions are 
evaluated as a shift of the data distribution with respect to reference distributions, and are 
measured with an unbinned likelihood fit calibrated with pseudo-data samples. Results 
from the two di-muon masses are combined in terms of a mean LEP CM energy shift 
giving 

Av^ = -334 ± 190(stat.) ± 76(syst.) MeV (1) 

where the main sources of systematic errors come from (i) possible biases in the muon 
polar angle measurement up to 88 = 0.2 mrad (52 MeV), (ii) the uncertainty from the 
shift fitting method (48 MeV) and (iii) from the absolute calibration of the muon momenta 
at the 0.5% level (24 MeV). This shift is 1.6cr from the nominal LEP CM average energy, 
consistent with no significant effect. 
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10 Combined results 



10.1 qqqq channel 

W mass 

For each of the standard, optimal CONE and PCUT reconstructions, the individual 
measurements of m w and T w at each CM energy are combined weighted by their 
statistical errors and systematic uncertainties as shown in Table 8. Correlations in these 
uncertainties with CM energy are taken into account. The W masses found from the one- 
parameter maximum likelihood fits to the data are given in Table 10. The experimental 

Table 10: W masses in the qqqq channel from all data for the standard, optimal PCUT and CONE 
reconstructions with corresponding statistical errors and systematic uncertainties (units are GeV/c 2 ). 



Reconstruction 


standard 


PCUT (3 GeV/c) 


CONE (R=0.4) 


Number of events 


4861 


4484 


4641 


m w 


80.481 


80.475 


80.502 


X 2 /dof 


7.4/7 


4.9/7 


5.1/7 


Statistical error 


0.058 


0.070 


0.070 


Experimental Uncertainty 


0.022 


0.028 


0.026 


FSI Uncertainty 


0.079 


0.028 


0.036 


Total error 


0.100 


0.081 


0.082 



systematic uncertainties are derived from all sources in quadrature including the LEP CM 
energy. The FSI uncertainties are the limits from colour reconnection and Bose-Einstein 
effects added in quadrature. Taking into account correlations, the W mass values are in 
good agreement for all three reconstructions. The total errors are also closely comparable 
but in the standard analysis, the systematic uncertainty due to FSI exceeds the statistical 
error. To suppress the dominant non-Gaussian contribution from colour reconnection, 
the result from the PCUT reconstruction with the smallest FSI uncertainty is selected to 
produce the final result to combine with the mass from the semileptonic channels. 

m§ = 80.475 ± 0.070(stat.) ± 0.028(syst.) ± 0.028(FSI) GeV/c 2 . 

The corresponding expected statistical error is 0.069 GeV/c 2 . 

The W masses with measured and expected statistical errors determined at each CM 
energy from the standard, CONE and PCUT analyses are given in Appendix B. 

W width 

The W widths found from the two-parameter maximum likelihood fits to the data are given 
in Table 1 1 for each reconstruction. The statistical error dominates the total systematic 
uncertainty in all three reconstructions. Therefore, the measurement from the standard 
analysis with the smallest total error is used in combination with the semileptonics to 
derive the most precise value for Tw- Thus, the W width from the qqqq channel is taken 
to be 

= 2.31 ± 0.12(stat.) ± 0.04(syst.) ± O.ll(FSI) GeV. 
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Table 11: W widths in the qqqq channel from all data for the standard, optimal PCUT and CONE 
reconstructions with corresponding statistical errors and systematic uncertainties (units in GeV). 



Reconstruction 


standard 


PCUT (3 GeV/c) 


CONE (R=0.4) 


r w 


2.31 


2.48 


2.34 


X 2 /dof 


9/7 


5/7 


16/7 


Statistical error 


0.12 


0.16 


0.15 


Experimental uncertainties 


0.04 


0.06 


0.05 


FSI uncertainties 


0.11 


0.03 


0.05 


Total error 


0.16 


0.17 


0.17 



The corresponding expected statistical error is 0.11 GeV. 

The W widths with measured and expected statistical errors determined at each CM 
energy from the standard analysis are given in Appendix C. 



10.2 ez/qq, /zz/qq and rz/qq channels 

W mass 



The mass values in the standard analysis from the one-parameter fits to the data, with 
the statistical and systematic errors including the LEP energy, are 

m ™w = 80.536 ± 0.087(stat.) ± 0.027(syst.) GeV/c 2 , 



rrv 



^qq 

w 



80.353 ± 0.082(stat.) ± 0.025(syst.) GeV/c 2 , 



m ™« = 80.394 ± 0.121(stat.) ±0.031(syst.) GeV/c 2 . 

The expected statistical errors are ±0.087, ±0.082 and ±0.122 GeV/c 2 for the e, fx and 
r semileptonic channels, respectively. 

The individual measurements of m-w for each channel are combined statistically at 
each CM energy. The combined semileptonic mw over all CM energies is determined by 
minimising a x 2 built from the full covariance matrix, taking into account all systematic 
uncertainties derived at each CM energy with the appropriate correlation and the 
statistical errors. The systematic uncertainties listed in Table 9 are found to be 89% 
correlated between ez/qq and /xz/qq, 85% between ez/qq and rz/qq and 89% between /zz/qq 
and rz/qq channels. 

The resulting combined mass for the semileptonic channels from the one-parameter 
fits is 

rrffii = 80.429 ± 0.054(stat.) ± 0.025(syst.) GeV/c 2 , 

with a x 2 /dof of 38/23. The expected statistical error is ± 0.054 GeV/c 2 . 

The W masses with measured and expected statistical errors determined at each CM 
energy are given in Appendix B. 

W width 

A two-parameter fit to the data gives the following results for the W width in the standard 
analysis for each channel: 

T^ qq = 1.84 ± 0.20(stat.) ±0.08(syst.) GeV, 
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T^ qq = 2.17±0.20(stat.) ± 0.06(syst.) GeV, 
r^ qq = 2.01 ± 0.32(stat.) ± 0.06(syst.) GeV, 

(2) 

where the expected errors are determined to be ±0.21, ±0.20 and ±0.31 GeV for the e, 
H and r channels respectively. The systematic uncertainties listed in Table 9 are found 
to be 100% correlated between ez/qq and /iz/qq, 43% between ez/qq and rz/qq and 48% 
between /zz/qq and rz/qq channels. 

The W widths with measured and expected statistical errors determined at each CM 
energy are given in Appendix C. 

The combined total width from the two-parameter fits in all £z/qq channels is 

Tw qq = 2.01 ± 0.13(stat.) ± 0.06(syst.) GeV, 
with a x 2 /dof of 15/21. The expected statistical error is ± 0.13 GeV. 

10.3 All channels 

The combined results from all channels using the PCUT results in the qqqq channel for 
the mass and standard results for the width are: 

m w = 80.444 ± 0.043(stat.) ± 0.024(syst.) ± 0.009(FSI) ± 0.009(LEP) GeV/c 2 , 
T w = 2.140 ± 0.090(stat.) ± 0.045(syst.) ± 0.046(FSI) ± 0.007(LEP) GeV. 

The combinations are performed in the same way as described in section 10.2. The x 2 /dof 
is 43/31 and 26/29 for the mass and width combinations, respectively. Alternatively, if 
the W mass from the CONE analysis in the qqqq channel is combined with those from 
the £z/qq channels, the mass is: 

m w = 80.453 ± 0.043(stat.) ± 0.023(syst.) ± O.Oll(FSI) ± 0.009(LEP) GeV/c 2 . 

Similarly, combining the W mass from the standard analysis in the qqqq channel with 
those from the £z/qq channels gives 

m w = 80.440 ± 0.043(stat.) ± 0.022(syst.) ± 0.019(FSI) ± 0.009(LEP) GeV/c 2 . 

To assess the effect of any unexpected correlation between the measured W mass and 
width on the one-parameter fits where the width is fixed to standard model values, the 
mass from each channel is compared with the corresponding two-parameter fit results. 
Combining all channels in the standard analysis, the difference is found to be 8 MeV/c 2 
indicating no significant effect. 

To investigate whether there is a significant difference between the masses from the 
qqqq and combined iuqq channels due to final state interactions, a fit is performed to 
extract this difference retaining all systematic uncertainties from Tables 8 and 9 except 
those from Bose-Einstein correlations and colour reconnection. The standard analysis in 
the qqqq channel is used to enhance any effects. The result is 

(mj™) - (tnj' qq ) = ±62 ± 76 (stat. ± syst.) MeV/c 2 , 
to be compared with the ±79 MeV/c 2 FSI uncertainty. 
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11 Conclusions 



The mass and width of the W boson have been measured from W pair events using 
the direct reconstruction of the invariant mass of their decay products in fully hadronic 
and semileptonic final states. Following constrained kinematic fits to each event, the W 
parameters were extracted by reweighting fully simulated invariant mass spectra to the 
measured distributions, employing an unbinned maximum likelihood procedure to find 
the best fits. To produce the most precise value of m^, one-parameter fits are performed 
where Tw varies with mw according to the Standard Model. Two-parameter fits, where 
Tow and Tw are allowed to vary independently, are used to measure the W width. 

All data collected at centre-of-mass energies between 183 and 209 GeV are fully 
reprocessed and analysed homogeneously to produce the final values with statistical errors. 
The systematic uncertainties are determined taking into account correlations between all 
channels and CM energies. For the W mass, these measurements are combined with the 
earlier published ALEPH results obtained from the total W pair cross sections at 161 [6] 
and 172 GeV [7] to produce the final result as follows: 

to w = 80.440 ± 0.043(stat.) ± 0.024(syst.) ± 0.009(FSI) ± 0.009(LEP) GeV/c 2 , 

where the first error is statistical, the second derived from all ALEPH experimental 
systematic uncertainties, the third from the final state Bose-Einstein and colour 
reconnection uncertainties in the qqqq channel and the last is the LEP energy uncertainty. 
The L3 and OPAL collaborations have recently published their results [45, 46] using all 
their available data. Also, earlier results have been published by DELPHI [47] as well as 
final results from the Tevatron Run I pp collider experiments using large samples of single 
W's decaying into electrons and muons [48]. 

No evidence is found for final state interactions between the W hadronic decay 
products in the qqqq channel. The limit on colour reconnection is derived from the 
search for any significant variation in the value of m w when low momentum particles or 
those between jets are progressively excluded in the invariant mass reconstructions. To 
minimise any colour reconnection effects, the W mass in the qqqq channel used in the 
final combination is taken from the reconstruction where all particles with momenta lower 
than 3 GeV/c are removed. 

This measurement of the W mass agrees with the earlier ALEPH measurement [5] 
and other measurements [45, 46, 47, 48] as well as with the indirect prediction from the 
Standard Model fit to electroweak observables [1]. The consistency with the electroweak 
fit is only possible if the Standard Model Higgs boson is light. 

Finally, from the 183-209 GeV data in all channels, the W width is determined to be 

T w = 2.14 ± 0.09(stat.) ± 0.04(syst.) ± 0.05(FSI) ± O.Ol(LEP) GeV, 
consistent with the other LEP measurements [45, 46, 47]. 
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Appendix A: Generator setup tunings used for CR 
studies 



Modified values of hadronisation and fragmentation parameters are tabulated for the 
model variants of GAL and ARIADNE with CR using the JETSET framework. The modified 
parameters with and without CR for the HERWIG model are also tabulated. The internal 
name of the parameters is given in each case. The unmodified parameters used in the 
generation of events for all models without CR are given in Ref. [49]. 





JETSET 


GAL 




ARIADNE 


Parameter 


standard and 




no CR 


intra-W only(AR21)/ 




SK models 




(AR20) 


intra and inter-W(AR2) 


azimuthal distribution 


3 





- 


- 


in PS MSTJ(46) 










momentum transverse 










width for hadron cr q t 


0.371 


0.364 


0.358 


0.352 


PARJ(21) (GeV) 










LUND fragmentation 










parameter b PAR J (42) 


0.805 


0.815 


0.823 


0.762 


Aqcd (GeV) PARJ(81) 


0.291 


0.307 






Q cut-off in PS (GeV) 
PARJ(82) 


1.52 


1.57 






strength parameter R 




0.039 






Colour Reconnection 










switch MSTA(35) 









1 (AR21) or 2 (AR2) 


Aqcd (GeV) PARA(l) 






0.230 


0.231 


ptmin cut (GeV) 
PARA (3) 






0.791 


0.781 


-Egiuon cut (GeV) 
PARA(28) 






0. 


2. 



Parameter 


HE 

no CR 


RWIG 

with CR 


Aqcd (GeV) QCDLAM 


0.190 


0.187 


Maximum cluster mass CLMAX 


3.39 


3.40 


Split cluster spectrum parameter 
light flavour clusters PSPLT(l) 
heavy flavour clusters PSPLT(2) 


0.945 
0.330 


0.886 
0.320 


Width of gaussian angle smearing CLSMR(l) 


0.58 


0.66 


Decuplet baryon weight DECWT 


0.71 


0.70 


CR probability PRECO 


0. 


1/9 


Gluon mass RMASS(13) 


0.774 


0.793 
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Appendix B: W masses from the standard, optimal 
CONE and PCUT analyses 



Table 12: Individual fitted tow values from the standard analysis for each channel and CM energy, 
together with the number of selected events and expected statistical errors. 



Channel 


CM energy 


Nevts 


mw 


Expected error 




(GeV) " 




(GeV/c 2 ) 


(GeV/c 2 ) 


qqqq 


183 


435 


80.525±0.168 


0.171 




189 


1169 


80.518±0.105 


0.102 




192 


234 


79.995±0.235 


0.235 




196 


556 


80.506±0.151 


0.151 




200 


627 


80.303±0.151 


0.146 




202 


283 


80.635±0.206 


0.218 




205 


589 


80.583±0.149 


0.146 




207 


968 


80.573±0.119 


0.118 


ez/qq 


183 


112 


80.440±0.265 


0.276 




189 


317 


80.437±0.170 


0.162 




192 


52 


80.621±0.502 


0.404 




196 


148 


80.420±0.244 


0.250 




200 


160 


80.607±0.247 


0.251 




202 


96 


80.203±0.303 


0.363 




205 


140 


81.089±0.276 


0.269 




207 


234 


80.620±0.218 


0.212 




183 


98 


79.991±0.265 


0.259 




189 


344 


80.185±0.160 


0.153 




192 


60 


80.483±0.385 


0.381 




196 


149 


81.109±0.246 


0.236 




200 


171 


79.884±0.237 


0.233 




202 


86 


81.210±0.324 


0.334 




205 


165 


80.409±0.250 


0.250 




207 


298 


80.277±0.186 


0.197 


rz/qq 


183 


97 


80.595±0.414 


0.396 




189 


306 


80.277±0.230 


0.232 




192 


59 


80.950±0.551 


0.548 




196 


158 


80.589±0.343 


0.347 




200 


163 


80.210±0.345 


0.348 




202 


70 


80.676±0.515 


0.509 




205 


149 


80.750±0.352 


0.370 




207 


224 


79.959±0.299 


0.296 
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Table 13: Individual fitted mw values from the CONE analysis in the 4q channel and CM energy, 
together with the number of selected events and expected statistical errors. 



Channel 


CM energy 
(GeV) " 


N ev t s 


m w 
(GeV/c 2 ) 


Expected error 
(GeV/c 2 ) 


qqqq 


183 


420 


80.607±0.214 


0.229 




189 


1116 


80.528±0.138 


0.132 




192 


224 


80.477±0.309 


0.319 




196 


537 


80.470±0.195 


0.194 




200 


589 


80.230±0.203 


0.191 




202 


272 


81.004±0.305 


0.278 




205 


565 


80.547±0.202 


0.197 




207 


918 


80.428±0.163 


0.159 



Table 14: Individual fitted mw values from the PCUT analysis in the 4q channel and CM energy 
together with the number of selected events and expected statistical errors. 



Channel 


CM energy 
(GeV) " 


N ev t s 


m w 
(GeV/c 2 ) 


Expected error 
(GeV/c 2 ) 


qqqq 


183 


409 


80.587±0.219 


0.230 




189 


1089 


80.529±0.138 


0.134 




192 


214 


79.935±0.316 


0.322 




196 


519 


80.517±0.197 


0.198 




200 


572 


80.357±0.203 


0.197 




202 


256 


80.614±0.320 


0.281 




205 


541 


80.333±0.205 


0.200 




207 


884 


80.588±0.165 


0.159 
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Appendix C: W widths from the standard analysis 



Table 15: Individual fitted T w values from two-parameter fits in the standard analysis for each channel 
and CM energy with expected statistical errors. The number of selected events is the same as the 
corresponding mass analysis in Appendix B. 



Channel 


CM energy 




Expected error 




(GeV) " 


(GeV) 


(GeV) 


qqqq 


X (JO 


X . (J iyT:_l_W . OOVJ 


yj .u(jw 




1 89 


2 ^20+0 242 

Zj . xj \J _1_ W . ZjtiZj 


214 

W . Zj It: 




192 


2 761 4-0 ^81 

Zj . 1 \J X _I_W .lJUl 


^30 




196 


2.022±0.331 


0.319 




200 


2 865±0 378 


0.311 




202 


9 1 34+n 545 


453 




205 

Zj VJ O 


1 790+0 345 


315 




207 


2 518+0 276 


247 


Ct^qq 


1 83 


9 4Q4+D 71 8 


n 6^7 




189 


1 662+0 353 


397 




1 92 

X u Zj 


Zj . X Ot:J L . X / Z, 


81 1 




196 


1.791+0.623 


0.611 




200 


2 909+0 715 


0.606 




202 


3 246+1 046 


764 




205 

Zj VJ O 


1 210+0 493 

X • Zj X \J _l W JO 


634 




207 

Zj VJ 1 


1 694+0 463 


508 


/ / 7/nn 


1 83 


1 894+0 630 

X . (J iyT:_l_W .UUU 


632 




1 89 


i 7054.(1 q^c. 

X. 1 iJO -i—\J .ooo 


fl 384 




192 

X <_/ Zj 




fl 783 




196 


2.667+0.700 


0.576 




200 


2.187+0.582 


0.572 




202 


2.566+0.957 


0.750 




205 


3.694+0.801 


0.606 




207 


2.186+0.454 


0.486 


ruqq 


183 


2.446+1.089 


0.848 




189 


1.720+0.527 


0.565 




192 




0.965 




196 


0.977+0.501 


0.776 




200 


1.808+0.761 


0.776 




202 


2.284+1.226 


0.925 




205 


2.171+0.830 


0.818 




207 


2.429+0.862 


0.700 



In the ruqq channel, the two-parameter fits fail to converge for the width at 192 and 
196 GeV even though the allowed range is 0.9 to 4.3 GeV. At 196 GeV, a one-parameter 
fit is successfully performed fixing the mass to the measured value given in Appendix B. 
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